Designing rigid carbon foams 
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We use ah initio density functional calculations to study the stability, elastic properties and elec- 
tronic structure of sp^ carbon minimal surfaces with negative Gaussian curvature, called schwarzites. 
We focus on two systems with cubic unit cells containing 152 and 200 carbon atoms, which are 
metallic and very rigid. The porous schwarzite structure allows for efficient and reversible doping 
by electron donors and acceptors, making it a promising candidate for the next generation of alkali 
ion batteries. We identify schwarzite structures that act as arrays of interconnected quantum spin 
dots or become magnetic when doped. We introduce two interpenetrating schwarzite structures that 
may find their use as the ultimate super-capacitor. 

PACS numbers: 81.05.Uw, 61.48.De, 73.22.-f, 71.20.-b, 62.25.-g 



I. INTRODUCTION 

Nanostructured carbon has become an intensely 
researched topic since the discovery of the Ceo 
moleculd^, a previously unknown carbon allotrope. 
Most research has focussed on structures with zero 
or positive Gaussian curvature such as fullerenes, 
nanotubes and graphene^^. Much less emphasis 
has been pla ced on structures with nega tive Gaus- 
sian curvaturj^l™ ^ l ^ l ^" l ^^ l ^^ l ^^ l ^^ l ^^ l ^^ l ^^ l ^^ l related to 
foams, which may be equally rigid and exhibit unusual 
electronic and magnetic properties. These systems, also 
called schwarzites, are space-filling contiguous structures 
formed of sp^ bonded carbon. They are unique in their 
ability to subdivide space into two disjoint, contiguous 
subspaces with labyrinthine morphology and the same 
infinite spatial extent. 

Here we focus on two schwarzite structures with cu- 
bic unit cells containing between 152 and 200 carbon 
atoms. We also postulate a structure with 400 atoms 
per unit cell as a representative of a class of interpene- 
trating, but disconnected schwarzite lattices. Using ah 
initio density functional calculations, we determine the 
equilibrium structure as well as the elastic, electronic and 
magnetic properties of pristine and electron or hole doped 
schwarzites. Our calculations indicate the possibility of 
tuning the concentration of K and CI atoms in order 
to achieve magnetic behavior. We find that the hole- 
doped C200 structure and similar systems may behave as 
an array of interconnected quantum spin dots. We dis- 
cuss potential energy storage applications of schwarzites, 
such as next-generation electrodes for alkali ion batter- 
ies, and calculate the capacitance of two interpenetrat- 
ing schwarzites, which may find their use as the ultimate 
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super-capacitor. 



II. COMPUTATIONAL METHOD 

Our ah initio calculations are based on tfie density 
functional tfieory (DFT) witfiin tfie local density approx- 
imation (LDA) and use tfie Ceperley-Alder excfiange- 
correlation functional^^ as parameterized by Perdew and 
Zungei^^. Interactions between valence electrons and 
ions are treated by norm-conserving pseudopotentialJ^ 
witfi separable non-local operators^^. Atomic orbitals 
witfi double-^ polarization are used to expand tfie elec- 
tron wave functions^^ witfi an energy cutoff of 210 Ry 
for tfie real-space mesfi. We use 0.02 Ry as tfie con- 
finement energy shift that defines the cutoff radii of the 
atomic orbitals. We sample the small Brillouin zone by 
8 /c-points in order to represent the Bloch wave functions 
for the momentum-space integration. All geometries are 
optimized using the conjugate gradient method^^, until 
none of the residual Hellmann-Feynman forces acting on 
any atom exceeds 1.56 x 10~^ Ry/ae, where is the 
Bohr radius. 



III. RESULTS 

A. Pristine schwarzite lattices 

Previous theoretical studies addressed schwarzites with 
24 (Ref. 9), 32 (Ref. 10), 48 (Ref.lIIJ, 56 fRefs. lTTI12t . 72 
(Ref.JJJ, 80 (Ref.JJ^, 168 (Refs. [HIISD, 192 (Refs.[ll 
[IS]), 200 (Ref.HH), 216 (Refs. HIl HH), and 224 (Ref. El) 
atoms per unit cell, or related general morphology is- 
sues (Refs. |T7l [18]). We focus on two structures with a 
primitive (P) minimal surface spanned by an underlying 
simple cubic lattice, which have not been studied previ- 
ously and which contain 152 and 200 carbon atoms per 
unit cell. The C152 structure is depicted in Fig. [TJa) 
and the C200 structure in Fig. [ij^b). The unit cells of 
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FIG. 1: Equilibrium geometry of P schwarzites with 152 (a) 
and 200 (b) carbon atoms per cubic unit cell with the lattice 
constant a and the volume V — . Total energy change 
as a function of V for the C152 (c) and C200 (d) schwarzites. 



FIG. 2: Electronic structure of schwarzites. Density of states 
of the optimized schwarzite structures with unit cells contain- 
ing 152 (a) and 200 (b) C atoms, smoothened by a Gaussian 
with the full width at half maximum of 0.03 eV. Correspond- 
ing band structure of the C152 (c) and C200 (d) schwarzites. 



either schwarzite contain one "core" structure with Oh 
symmetry and a negative local Gaussian curvature. The 
six extremities of the cores are connected to neighbor- 
ing cores by (4, 4) carbon nanotube segments of various 
length to form a contiguous lattice. Negative Gaussian 
curvature is introduced by the presence of heptagons in 
the graphitic lattice. Our systems contain exactly 24 
heptagons per unit cell in accordance with Euler's the- 
orem that, along with 24 hexagons, form the cores of 
the schwarzite. Each interconnect contains one nanotube 
unit cell with 8 hexagons in the C152 structure and two 
unit cells with 16 hexagons in the C200 structure. 

To identify the equilibrium structure and elastic prop- 
erties of the two schwarzite structures, we performed a 
series of structure optimizations at fixed unit cell volume. 
The total energy change /S.E{y) with respect to the op- 
timum structure is presented for the two structures in 
Figs. [TJc) and [ijd) along with Murnaghan function fits 
of our data. For the C152 schwarzite we find the binding 
energy Ecoh = 10.50 eV/atom at the equilibrium vol- 
ume Veq = 2515 A^, corresponding to the lattice con- 
stant a = 13.60 A, gravimetric density p = 1.21 g/cm^ 
and bulk modulus B = 111 GPa. The less dense 
C200 schwarzite structure is slightly more stable with 
Ecoh = 10.52 eV/atom at Veq = 4141 A^, correspond- 
ing to a = 16.06 A, but has a smaller gravimetric den- 
sity p = 0.96 g/cm^ and bulk modulus B = 76 A GPa. 
The reduction of the binding energy with respect to the 
value of 10.88 eV/atom in a graphene monolayer is caused 
by the strain associated with curvature and is similar in 
value to carbon nanotubes. We find both structures to 
be substantially lighter than graphite {p = 2.27 g/cm^) 
and diamond^^ {p = 3.54 g/cm^). Also the bulk modulus 
of the stiffer C152 is about four times smaller than the 



diamond valu^ Bexpt = 442 GPa. 

To characterize the electronic properties of these 
schwarzites, we calculated their electronic density of 
states and band structure. The electronic densities of 
states of Ci52 and C200, presented in Figs.|2|a) and|2|b), 
indicate that both structures are metallic. The origin of 
the individual peaks in the density of states can be traced 
back to the band structure graphs in Figs.[2|c) and|2|d). 
The schwarzites exhibit narrow band gaps above and, in 
the case of C2005 ^Iso below the Fermi level. Further- 
more, as seen in Figs. |2|b) and[2jd), the C200 structure 
exhibits a 3-fold degenerate flat band about 0.62 eV be- 
low the Fermi level. Since in neither structure the Fermi 
level falls into a sharp peak, both structures are non- 
magnetic. 

Besides the nonzero value of the density of states at Ep 
in the schwarzites we study, which makes them metal- 
lic, also the nature of those states is very important for 
the dielectric response of the system. Therefore, we plot 
the electron density associated with the frontier states at 
the Fermi level of C152 in Fig. jsja) and that of C200 in 
Fig. [sj^b). We find the frontier states to be highly delo- 
calized in both systems. In the larger C200 schwarzite, 
we note that the frontier states have a node in the middle 
of the tubular interconnects between the cores, indicated 
by an arrow in Fig. [3|b). As seen in Fig. [Sjc), which 
depicts the electron density of C152 associated with the 
~0.8 eV wide band mostly below the Fermi level, the cor- 
responding states are distributed nearly uniformly across 
the structure. To understand the origin of the flat band 
in the C200 schwarzite 0.62 eV below the Fermi level, we 
present the associated electron density in Fig. [sj^d). We 
may characterize this state as an array of non-interacting 




FIG. 3: (Color online) Electron densities in schwarzite struc- 
tures containing 152 and 200 atoms per unit cell, presented 
as isodensity surfaces. Electron density associated with fron- 
tier states of Ci52 (a) and C200 (b) in the energy range 
Ef - 0.15 eV< E < Ef + 0.15 eV, presented at the iso- 
density value oflO~^A~^. (c) Electron density of C152 as- 
sociated with states in the energy range Ef — 0.72 eV< E < 
Ef + 0.10 eV, presented at the isodensity value of 10~^ A~^. 
(d) Electron density of C200 associated with the localized 
state near Ef — 0.62 eV, presented at the isodensity value 

of 10-^ A-^ 

quantum dots centered in the negative Gaussian curva- 
ture regions at the cores of the C200 structure. We expect 
such a state to occur also in other schwarzite structures of 
the same type, which have larger unit cells due to longer 
tubular interconnects. 



B. Interpenetrating schwarzite lattices 

The open space in the C152 and C200 structures, de- 
picted in Figs.jlja) andjljb), can be filled by other atoms 
as external dopants. As we show in the following and de- 
pict in Figs. Qa) and[4]^b), the open space in the C200 
schwarzite is wide enough and ideally suited to accom- 
modate even an identical replica of the C200 schwarzite 
in an unstrained lattice. 

In either schwarzite described in this manuscript, the 
tubular interconnects are composed of (4, 4) carbon nan- 
otubes. As seen in Figs. [4]^a) and|4]^b), the closest ap- 
proach distance between the sublattices occurs in the 
geometry depicted in the inset of Fig. |4]^c), with two 
(4, 4) nanotube segments standing normal to each other. 
To estimate the optimum inter-wall separation, we cal- 
culated the inter-tube interaction in this geometry and 
plotted the data points along with a Morse function fit 
in Fig. |4]^c). We found the optimum inter-tube distance 
to be deq = 7.8 A, corresponding to the closest inter- 
wall distance of 2.4 A, somewhat smaller than the inter- 
layer distance in graphite. The corresponding distance 
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FIG. 4: (Color online) Equilibrium structure of two interpen- 
etrating C200 schwarzite lattices, shown in Fig.[T]^b), with 400 
C atoms in total per unit cell, shown in perspective (a) and 
top (b) view, (c) Interaction energy between two (4, 4) nan- 
otubes, with the geometry shown in the inset, representing 
the interaction between the tubular interconnects in the two 
interpenetrating C200 schwarzites. (d) Electronic density of 
states of the C400 schwarzite. 



between adjacent tubular interconnects of the C200 sub- 
structures in C400, as depicted in Figs. [4]^ a) and|4jb), is 
d = 8.00 A, which indeed is very close to the optimum 
inter-tube distance value deq = 7.8 A. Since the C400 
schwarzite is structurally a superposition of two C200 
schwarzite lattices, its gravimetric density is twice that 
of C2005 namely p = 1.92 g/cm^. Also the bulk modulus 
of C400, B = 153 GPa, is twice that C2007 which is about 
one third of the diamond value. 

As the electronic density of states of the C200 sublat- 
tices, also that of the C400 schwarzite, shown in Fig.|4]^d), 
has a large density of states at the Fermi level, indicat- 
ing a metallic system. Similar to C2005 also C400 has a 
narrow band about 0.6 eV below Ep^ associated with a 
lattice of very weakly coupled quantum dots. Careful 
comparison of Figs.[2|b) and[4]^d) indicates that the cen- 
tral band around Ep becomes wider in the C400 lattice in 
comparison to the C200 lattice. This band broadening is 
accompanied by a corresponding narrowing of the band 
gaps above and below this central band. 



C. Intercalated schwarzites 

As mentioned above, the empty space delimited by the 
schwarzite minimal surface can accommodate intercalant 
ions in order to shift the Fermi level. We consider K as 
a model donor that has been used widely in graphite 
intercalations compounds (GICs)P3. Due to its highest 
electronic affinity in the periodic table^^ of 3.61 eV, we 
choose CI as a model acceptor. 



As in the previous subsection, we focus on the C200 
schwarzite structure. To maximize charge transfer be- 
tween the inter calant atoms and the schwarzite lattice, 
we place one K or CI atom in the center of each core of 
the schwarzite. To estimate the charge transfer between 
the intercalants and the schwarzite, we performed a Mul- 
liken population analysis and, alternatively, integrated 
the total electron density within a sphere surrounding 
the inter calant ion. We found the Mulliken charge on K 
to be +0.96e and that on CI to be — 0.42e. Electron den- 
sity integration around the intercalant within schwarzite 
only makes sense up to the size of the cavity inside the 
core structure. Considering a sphere with a radius of 
R = 2A A around CI yields 7.00 valence electrons in case 
of an isolated atom and 7.42 valence electrons in case of 
Cl@C2oo- Since the atomic radius of K is much larger, an 
isolated atom contains only 0.56 valence electrons within 
a sphere with R = 2.4 A. Upon intercalation, in K@C2oo, 
it loses almost all of the valence charge, so that only 0.01 
valence electrons remain within this sphere. 

By subtracting the total energy of C200 and the iso- 
lated atom A from that of A@C200 7 where A=K or CI, we 
found that the K atom stabilizes the unit cell by 2.080 eV, 
whereas the CI atom provides a smaller stabilization by 
0.636 eV per unit cell, partly due to the smaller charge 
transfer. 

The ability to efficiently transfer most of the valence 
charge from the alkali atom to the graphitic structure 
is used currently in Li ion batteries, where Li is in- 
tercalated in-between graphene layers in graphite, caus- 
ing an ^30% expansion in the direction normal to the 
graphene layers. Repeating expansion and contraction 
during charge/discharge cycles is known to reduce the 
lifetime of alkali ion batteries, as it causes structural 
rearrangement of the grains and eventually blocks dif- 
fusion pathways. No such adverse effects should occur 
in the rigid, but still conductive schwarzite structures, 
which may yield a new generation of advanced alkali ion 
batteries. As a different potential application of doped 
schwarzites, we wish to mention their potential for hydro- 
gen storage, similar to Ca doped Ceo crystals that proved 
to store reversibly molecular hydrogerP^up to 8.4% (wt). 

As suggested by the smaller amount of charge trans- 
ferred from the C200 schwarzite to CI ions, the effect of 
acceptor intercalation is much smaller than that of donor 
intercalation. Nevertheless, hole doping should make it 
possible to shift Ep down. Especially interesting appears 
the possibility to move the Fermi level into the region of 
the threefold degenerate flat band. Partial occupation 
of the band is expected to change the system to a fer- 
romagnet with a magnetic moment per unit cell up to 
6/j.B' To determine the necessary amount of doping, we 
estimated the occupied portion of the band near Ep to 
carry 6 electrons. Combining this with the fact that the 
flat band also contains 6 electrons, we conclude that hole 
doping C200 with 7—11 electrons per unit cell should, 
assuming that the rigid band model applies, move the 
Fermi level into the region of the flat band. Using the 
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FIG. 5: (Color online) Electron density difference 
An(A@C2oo) = ntot{A<^C2oo) - utotiA) - ntot{C2oo), pro- 
viding insight into the charge flow in the C200 schwarzite fol- 
lowing intercalation with A = K, CI atoms, (a) An{K(^C2oo) 
at the isodensity value +0.04 A~^. (b) An(C/@C2oo) at the 
isodensity value —0.04 A~^. 



charge transfer estimated for one CI atom in C2007 this 
could be achieved by intercalating between 10 — 30 CI 
atoms per unit cell. Using interatomic distances from 
free CI2, we believe that in the optimum case, Cle clus- 
ters could be accommodated in the cores. Assuming that 
also each tubular interconnect could accommodate up to 
two CI atoms, the CI doped system should barely reach 
the criterion of partly depleting the initially flat band by 
CI doping, causing ferromagnetism to occur. 

To get an impression of how uniformly the transferred 
charge is distributed across the schwarzite lattice, we plot 
in Fig.jsja) the distribution of the excess electron density 
in K@C2oo and in Fig. [5|b) the corresponding distribu- 
tion of the deficient electron density in Cl@C2oo- Since 
both ions are well separated from the schwarzite lattice, 
we find that the transferred charge is distributed rather 
uniformly across the schwarzite in both cases. 

Also the two interpenetrating C200 lattices, forming 
the C400 schwarzites, can be differentially doped, one 
with K and the other with CI atoms, giving a net posi- 
tive charge to one and a net negative charge to the other 
sublattice, similar to a capacitor. We superposed the 
K@C2oo and Cl@C2oo lattices and found a net stabi- 
lization with respect to separated K@C2oo and Cl@C2oo 
doped schwarzites of —15.735 eV. Similar to the doped 
sublattices, the Mulliken population analysis indicates 
that the net charge on K is close to +0.70e and that 
on CI is close to — 0.42e. The opposite charges are then 
accommodated in the adjacent C200 substructure. 

In the following we will estimate the capacitance of 
the schwarzite structure per 1 cm^ volume. The capac- 
itance of a parallel-plate capacitor of area A and inter- 
plate distance (i, defining the volume V = Ad^ is given 
by C = eoA/d in SI units. For the sake of reference, 
an idealized capacitor with V = 1 cm^, formed by a 
graphene bilayer with d = 0.335 nm, would have the 
area A = 3.0x10^ m^ and a capacitance of 79.3 F. 
In the unit cell of the C400 schwarzite with the vol- 
ume V = 4.141 nm^, 200 carbon atoms are associated 
with one electrode. With an estimated area per carbon 
atom of 2.62 A^, same as in graphene, we estimate that 
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the area of one electrode per cm^ schwarzite should be 
1.3x10^ m^. Assuming furthermore the same average in- 
terlayer distance as in a graphene bilayer, we estimate 
the capacitance of the C400 schwarzite to be 34.3 F per 
cm^ material. 

As an independent way to estimate the capacitance 
of the C400 schwarzite, we carefully inspected the total 
energy differences between neutral C200 ^ind C400 sys- 
tems, as well as the total energies of K@C200 7 Cl@C2oo 
and C400 containing both K and CI. The latter struc- 
ture, which contains one positively and one negatively 
charged C200 minimal surface per unit cell, is a superpo- 
sition of K@C2oo and Cl@C2oo- The stabilizing interac- 
tion between the charged C200 substructures of 15.735 eV 
is partly due to the chemical interaction between the 
two electrodes, which we estimated to be 15.464 eV per 
unit cell in the pristine system. We assign the remain- 
ing part of the interaction energy, 0.271 eV, to be the 
energy associated with the electric field between the ca- 
pacitor plates. To understand the meaning of this energy 
gain, we consider two parallel plate electrodes carrying 
constant charges -\-Q and — Q, approaching from infin- 
ity to form a parallel-plate capacitor of capacitance C, 
thereby gaining the energy AU = Q'^/{2C). Consider- 
ing AU = 0.271 eV and \Q\ = 0.42e, the smaller of the 
charges transferred between the intercalant atoms and 
the neighboring electrodes, we can estimate the capac- 
itance C per unit cell and the corresponding value in 
the bulk C400 schwarzite material, which turns out to be 
12.6 F per cm^. 

We would like to point out that the above capaci- 
tance values are based on rough approximations and thus 
should be considered as order-of-magnitude estimates. 
Still, we are pleased that our estimate based on total 
energy differences in the intercalated system lies close to 
the estimate based on area and inter-plate distance. In 



any case, the estimated capacitance exceeds that of nowa- 
days' capacitors by several orders of magnitude. Should 
a material close to the postulated C400 schwarzite ever be 
synthesized, we must also consider the need to connect 
the two interpenetrating sublattices to two leads, which 
will be a nontrivial task. We also need to point out that 
the close proximity of the two electrodes in the schwarzite 
material promotes tunneling, which would eventually dis- 
charge the electrodes over time. Still, the expected ca- 
pacitance of tens of Farads per cubic centimeter is an 
appealing prospect. 

In summary, we studied the stability, elastic proper- 
ties and electronic structure of sp'^ carbon periodic min- 
imal surfaces with negative Gaussian curvature, called 
schwarzites, using ab initio density functional calcula- 
tions. Our studies of two primitive minimal surfaces, 
spanned by an underlying simple cubic lattice with 152 
and 200 carbon atoms, indicate that these systems are 
very stable, rigid, and electrically conductive. The 
porous schwarzite structure allows for efficient and re- 
versible doping by intercalation of electron donors or ac- 
ceptors. We identified systems that should act as arrays 
of interconnected quantum spin dots or, when doped, 
exhibit ferromagnetic behavior. We also introduce two 
interpenetrating schwarzite structures as an unusual sys- 
tem that may find its use as the ultimate super-capacitor. 
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ported by NRF of Korea grant KRF-2009-0074951. 
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